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ABSTRACT 

Context. Intermediate-mass young stellar objects (YSOs) provide a link to understand how feedback from shocks and UV radiation 
scales from low to high-mass star forming regions. 

Aims. Our aim is to analyze excitation of CO and H 2 0 in deeply-embedded intermediate-mass YSOs and compare with low-mass 
and high-mass YSOs. 

Methods. Herschel/PACS spectral maps are analyzed for 6 YSOs with bolometric luminosities of Lb 0 i ~ 10 2 - 10 3 L 0 . The maps 
cover spatial scales of ~ 10 4 AU in several CO and H 2 0 lines located in the ~ 55 - 210 pm range. 

Results. Rotational diagrams of CO show two temperature components at T mt ~ 320 K and 7j- ot ~ 700 - 800 K, comparable to low- 
and high-mass protostars probed at similar spatial scales. The diagrams for H 2 0 show a single component at T'rot ~ 130 K, as seen 
in low-mass protostars, and about 100 K lower than in high-mass protostars. Since the uncertainties in T Iot are of the same order as 
the difference between the intermediate and high-mass protostars, we cannot conclude whether the change in rotational temperature 
occurs at a specific luminosity, or whether the change is more gradual from low- to high-mass YSOs. 

Conclusions. Molecular excitation in intermediate-mass protostars is comparable to the central 10 3 AU of low-mass protostars and 
consistent within the uncertainties with the high-mass protostars probed at 3 ■ 10 3 AU scales, suggesting similar shock conditions in 
all those sources. 
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1. Introduction 


Feedback processes associated with the collapse of protostel- 
lar envelopes at 10 3 - 10 4 AU scales limit the accretion onto 
the protostar and contribute to the overall low efficiency of 
transferring gas into stars on global scales (lOffner et alJl2009t 
iKrumholz et alJ[2014l) . Calculating the physical conditions help 
to identify the most relevant phenomena and constrain their role 
in the star formation process dEvansll 19991) . Since gas in proto- 
stellar envelopes is heated to temperatures much higher than the 
dust temperatures, molecular transitions are the suitable tracers 
of physical conditions of hot (T > 100 K) gas around proto¬ 
stars. In particular, the far-infrared (IR) lines of CO and H 2 0 
dominate the cooling of hot and dense gas (iGoldsmith & Langeii 
119781 ). The excitation of CO and HiO depends on the local phys¬ 
ical conditions (temperature, density) and thus is crucial to de¬ 
termine which physical mechanisms are responsible for the gas 
heating and to study whether the energetics involved in the feed¬ 
back scale from low- to intermediate- to high-mass young stellar 
objects (YSOs). 


Recent observations of CO and H 2 0 lines with the 


Photodetector Arra 


I 


Camera and Spectrometer (PACS, 
on board Herschel found evidence for 
10 4 cm‘ 3 ) and hot (> 300 K) gas 


iPoglitsch et al.l 1201 
large columns of dense (> 

towar d s low-mass (Lh n i 5, 1 0 2 L o ) protosta r s (Iva n Kempen et al 


1201 ft iHerczeg et all 120121 : iManoi et al.1 12013 : Karska et al 


120131: iGreen et alll2013t iLindberg et al.ll2014h . originate largely 
from UV-irradiated shocks associated with jets and winds 
( Karska et al. 2014bl) . CO and H 2 0 line luminosities of the 
high-mass protostars (with Lboi ~ 10 4 — 10 6 L 0 ) follow 

the correlations with bolometric luminosities found in the 
low-mass protostars (iKarska et al.l l2014al) and show similar 
velocity-re solved line p r ofiles r egardless of the mass of the 
protostar jYildi zeLaL 120131: ISan Jose-Garcla et al.l 12013b 
Ivan der Tak et al.l l2013h . In contrast, rotational temperatures of 
H 2 0 are lower and the H 2 0 fraction contributed to the total 
cooling in lines with respect to CO is higher for the low -mass 
protostars (IKarska et al.l 1201 4at iGoicoechea et al.l 1201 5l) . sug¬ 
gesting that the physical mechanism causing the excitation in 
low- and high-mass protostars are different. 

Intermediate-mass YSOs (withLboi ~ 10 2 - 10 3 L .|]]) provide 
a natural link between low- and high mass protostars, but their 
far-IR CO and H 2 0 emission has only been studied for a single 
protostar position of NGC 7129 FIRS 2 dF ich et al. l2010l) and the 
outflow position of NGC 2071 (iNeufeld et al. 20141) . CO emis¬ 
sion alone was analyzed for two intermediate-mass protostars 


1 We use bolometric luminosity as a proxy of the protostellar mass 

for practical and historical reasons, but we note that Lboi changes sig- 

nificantly during the protostellar phase if the accretion is episodic 

lYoun g & Evans 2005; Dunham et al. 2010). Moreover, some of our 

intermediate-mass sources may be in fact a collection of unresolved 

low-mass protostars. 


1 


































































M. Matuszak et al. 2015: Intermediate-mass YSOs with PACS (RN) 





NGC2071 J 

310(30) K 



CO 700(50) 


■ I - - - - I ■ ■■■ I - 


v 


NGC7129 : 

370(45) K 



CO o 

710(70) K 


■ i - - - - i ■ ■■■ i - 


2000 4000 


Ejk B (K) 




Fig. 1. Rotational diagrams of CO. The base-10 
logarithm of the number of emitting molecules 
from the upper level, Af up , divided by the de¬ 
generacy of the level, g up , is shown as a func¬ 
tion of energy of the upper level in kelvins, E up . 
Detections are shown as filled circles, whereas 
three-sigma upper limits are shown as empty 
circles. Empty upper triangle corresponds to 
the line flux calculated using a smaller area on 
the map than the rest of the lines. Blue lines 
show linear fits to the data and the correspond¬ 
ing rotational temperatures. Errors associated 
with the least-square linear fit are shown in 
brackets. 


Fig. 2. Similar to Figure[T|but for TEO. 


in Orion, HOPS 288 and 370 dManoi et al.11201 3 ). In this paper, 
we present the analysis of PACS spectra for the full sample of 
intermediate-mass protostars from the ‘Water in star forming re¬ 
gions with Herschel’ (WISH) key program (van Dis hoeck et alJ 
1201 ill , including the maps of NGC 7129 and NGC 2071 centered 
on the YSO position. These results complement the work by 
IWamofler et al .i ( 201 3). which describes the OH excitation in our 
source sample and the sample of low- and high-m ass p rotostars 
for which CO and H 2 O emission is discussed in Karsk a et al.l 
( 2013 1 2014a ). The main question addressed is whether CO and 
H 2 O rotational temperatures differ from low- to high-mass pro¬ 
tostars. 

The paper is organized as follows. §2 briefly introduces the 
observations, §3 the excitation analysis using rotational dia¬ 
grams, and §4 discusses the results. 


2. Observations 

Our sample includes 6 YSOs with bolometric luminosities from 
70 to 2000 L q and located at an average distance of 700 pc (see 
Table |T}. The sources were selected based on their small dis¬ 
tances (< 1 kpc) and location accessible for follow-up observa- 
ti ons from the southern hem isphere (for more details see §4.4.2. 
in Ivan Dishoeck et alJi201 lh . Spectroscopy for all sources was 
obtained with PACS as part of the WISH program. For observ¬ 
ing details see Table lA.fl in the Appendix. 

With PACS, we obtained single footprint spectral maps cov¬ 
ering a field of view of ~ 47" x 47" and resolved into 5x5 array 
of spatial pixels (spaxels) of ~ 9.4" x 9.4" each. At the dis¬ 
tance to the sources, the full array corresponds to spatial scales 
of~2-6xl0 4 AU in diameter, of the order of full maps of low- 
mass YSOs (~ 10 4 AU) in iKarska et al.l d2013h and the central 
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Table 1 . CO and HiO rotational excitation 


Source 

D 

(pc) 

^bol 

0®) 

Warm CO 
r rot (K) log W N 

Hot CO 

r r „,(K) log ln Af 

HiO 

Trot(K) log,„)V 

AFGL 490 

1000 

2000 

300(20) 

51.6(0.1) 



90(100) 

48.3(1.2) 

L 1641 S3MMS1 

465 

70 

325(30) 

49.9(0.1) 



150(100) 

46.8(0.1) 

NGC 2071 

422 

520 

310(30) 

51.2(0.1) 

700(50) 

50.1(0.1) 

135(100) 

47.9(0.1) 

Vela 17 

700 

715 

265(25) 

50.8(0.2) 



90(40) 

47.8(0.5) 

Vela 19 

700 

776 

335(50) 

50.4(0.2) 



130(80) 

47.2(0.4) 

NGC 7129 FIRS 2 

1250 

430 

370(45) 

50.7(0.2) 

710(70) 

49.9(0.2) 

130(50) 

47.8(0.1) 


Notes. Distances and bolometric luminosities are taken from lWampfler et aD <12013l) and references therein. 


spaxel spectra of more distant high-mass YSOs (~ 3 x 10 4 AU) 
in lKarska et al.l (l2014al) . 

The observations were taken in the line spectroscopy mode, 
which provide deep integrations of 0.5-2 gm wide spectral re¬ 
gions within the ~55-210 pm PACS range. Two nod positions 
were used for chopping 3' on each side of the source (for de¬ 
tails of our observing strategy and basic reduction methods, see 
iKarska et al.112013l) . The full list of targeted CO and PDO lines 
and the calculated line fluxes are shown in Table IA.2I in the 
Appendix. The quality of the spectra is illustrated in Figure IaTI 
in the Appendix. We note that the simultaneous non-detections 
of the H 2 O 7 16 -7ov line at 84.7 fin 1 and detections of the FFO 
8i8-7 07 line at 63.3 /um in L 1641, NGC 2071, and NGC 7129 
FIRS 2 are related to the structure of the HiO energy levels and 
not the differences in the sensitivities of the instrument in the 
2nd and 3rd order observations. The 8 is-7o 7 line is a backbone 
transition, with the Einstein A coefficient higher than the one for 
the 7i6-7 07 line. 

The data reduction was done using HIPE v.13 with 
Calibration Tree 65 and subsequent analysis with customized 
IDL programs (see e.g. iKarska et alJl2014b ). The fluxes were 
calculated using the emission from the entire maps. Figure fACl in 
the Appendix illustrates that both the line and continuum emis¬ 
sion peaks approximately at the source position, with small shifts 
in continuum due to mispointing. The extent of the line emission 
follows typically the continuum pattern with the exception of 
Vela 17 where the line emission extends from NE to SW direc¬ 
tion, while the continuum is centrally peaked. There is no con¬ 
tamination detected from the nearby sources or their outflows in 
the targeted regions. 



Lboi (L©) 

Fig. 3. Rotational temperatures of ‘warm’ (£,’ up = 14 - 

24) CO (top) and FFO (bottom) as function of bolomet- 
ric luminosity. Orange cir cles show YSOs from jGreen et al.1 
(l2013h . IKarska et ail n2()\ 3l) . and lManoi et alJ (120131) . Dark blue 
crosses show intermediate-mass YSOs from the W ISH program 
and light blue diamonds high-mass YSOs from IKarska et al.l 
( 2014a’) . A few sources with luminosities ~ 10 2 L (:) shown in 
orange crosses are from iManoi et al.l d2013i) and can be also re¬ 
garded as intermediate-mass YSOs. Uncertainties of the CO and 
H 2 O rotational temperatures are shown for L1641 S3 MMS1, 
which is representative of the sample. 


3. Rotational diagrams 

3.1. Results 

Figure Q] shows rotational dia grams of CO calculated for all 
sources in the same way as in lKarska et al.l (2 014 a). The cor¬ 
responding rotational temperatures, T mt , and total numbers of 
emitting molecules, N, are shown in Table [T] 

All sources show a 300 K, ‘warm’ CO comp onent 
(IManoi et al.ll2013l IKarska et al.ll2013l iGreen et aDl2013l) . with 
a mean temperature of ~ 320 K for the range of total numbers 
of emitting molecules, N ~ 10 50 - 10 51 . In addition, NGC2071 
and NGC7129 show a ‘hot’ CO component with temperatures of 
~ 700 K and N ~ 10 50 . 

The rotational diagrams for H 2 O, presented in Figure [2] 
show a single component with a mean temperature of ~ 120 
K. The scatter due to subthermal excitation and high opacities 
exceeds the uncertainties in the observed fluxes, similar to dia¬ 
grams of low- and high-mass YSOs (for the discussion of both 
effects see §4.2.2 of IKarska et alJ l2014a l. The corresponding 


numbers of emitting molecules are about 3 orders of magnitude 
lower than for CO, N ~ 10 47 - 10 48 . 


3.2. Comparison to low- and high-mass sources 

Figure [3] shows a comparison between rotational temperatures 
obtained for CO and H 2 O for the intermediate-mass sources pre¬ 
sented here and these quantities determined in the same way for 
low- and high-mass YSOs. The comparison is restricted only to 
the ‘warm’ component seen in CO rotational diagrams due to the 
low number of sources with detections of the ‘hot’ CO compo¬ 
nent; those will be discussed in detail in Karska et al. (in prep.). 
A bolometric luminosity is used here as a proxy for the proto- 
stellar mass, but in fact some of our intermediate-mass sources 
may be a collection of unresolved low-mass protostars. 

The median T rot of CO in low- mass pr otostars is 325 K, using 
the results from the WISH ( Karska et alJl2013l) . ‘Dust, Ice, and 
Gas in Time’ dGreen et alJl2013 ). and ‘Herschel Orion Protostar 
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Survey’ dManoi et al.l 12013 1) programs for a total of about 50 
sources. A comparable value of ~290 K was found for 10 high- 
mass sources (Lboi ~ 10 4 - 10 6 L e ) in iKarska etaD (l2014ai) . 
The range of CO T rot of 265-370 K (see Table Q} determined 
here for intermediate-mass YSOs is thus fully consistent with 
previous results. The fact that the 300 K CO component does 
not depend on the source bolometric luminosity over 6 orders of 
magnitude suggests the origin in a shock associated with the jet / 
winds impact on the envelope rather than in a photodissociation 
region where IisL s hould scale with the UV flux and luminosity 
( vanKemgeneUd. 2010t Visseret al.|2012|: Manoi et a j1l2013t 
iFlower & Pineau des Foretsll2013t kristensen et al.ll20i*3l) . 


The median FFO rotational temperatures (see the bottom 
panel of Figure [3} for low- a nd high-mass YSOs are ~ 130 
K and ~ 230 K, respectively (IKarska et al.l 1201 3l 1201 4al) . The 
range of temperatures obtained for the intermediate-mass YSOs 
(90-150 K, 'Fable |Tj is thus comparable to the low-mass sources. 
However, the uncertainties in the rotational temperatures are 
high, of the order of ~ 100 K, and do not account for the opti¬ 
cal depth effects, the density effects, and the possible complexity 
of the line profiles (absorption and emission components). It is 
therefore unclear if there is a true jump in rotational tempera¬ 
tures at Lboi ~ 10 4 L q , or a smooth trend toward higher values 
of 7j ot . In either case, these higher excitation temperatures could 
be due to higher densities in the more massive envelopes. 


Karska, A., Herpin, F., Bruderer, S., et al. 2014a, A&A, 562, A45 
Karska, A., Kristensen, L. E., van Dishoeck, E. F., et al. 2014b, A&A, 572, A9 
Kristensen, L. E., van Dishoeck, E. F.. Benz, A. O., et al. 2013, A&A, 557, A23 
Krumholz, M. R.. Bate, M. R., Arce, FI. G., et al. 2014, Protostars and Planets 
VI, University of Arizona Press (2014), eds. FI. Beuther, R. Klessen, C. 
Dullemond, Th. Flenning 

Lindberg. J. E., J0rgensen, J. K„ Green, J. D., et al. 2014, A&A, 565, A29 
Manoj, P, Watson, D. M., Neufeld, D. A., et al. 2013, ApJ, 763, 83 
Mottram, J. C., Kristensen, L. E., van Dishoeck, E. F., et al. 2014, A&A, 572, 
A21 

Muller, FI. S. P, Schloder, F.. Stutzki, J., & Winnewisser, G. 2005, Journal of 
Molecular Structure, 742, 215 

Muller, H. S. P, Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, 
L49 

Neufeld, D. A., Gusdorf, A.. Glisten. R., et al. 2014, ApJ, 781, 102 
Offner, S. S. R., Klein, R. I., McKee, C. F., & Krumholz, M. R. 2009, ApJ, 703, 
131 

Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, 
J. Quant. Spec. Radiat. Transf., 60, 883 
Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 
San Jose-Garcla, I., Mottram, J. C., Kristensen, L. E., et al. 2013, A&A, 553, 
A125 

van der Tak, F. F. S., Chavarria, L., Herpin, F., et al. 2013, A&A, 554, A83 
van Dishoeck, E. F., Kristensen, L. E., Benz, A. O., et al. 2011. PASP, 123, 138 
van Kempen, T. A., Kristensen, L. E., Herczeg, G. J., et al. 2010, A&A, 518, 
L121 

Visser, R., Kristensen, L. E., Bruderer, S., et al. 2012, A&A, 537. A55 
Wampfler, S. F., Bruderer, S., Karska, A., et al. 2013. A&A, 552, A56 
Ytldtz, U. A.. Kristensen, L. E., van Dishoeck, E. F., et al. 2013, A&A, 556, A89 
Young, C. H. & Evans. II, N. J. 2005, ApJ. 627, 293 


4. Summary 

We analyze the excitation of far-infrared CO and H 2 O lines in 
6 intermediate-mass YSOs observed in the WISH survey and 
compare the results to low- and high-mass protostars. Rotational 
temperatures of CO and H 2 O are found to be ~ 320 K and 
~ 120 K, respectively, and are consistent with low-mass and 
high-mass YSOs within the uncertainties. The large uncertain¬ 
ties in the H 2 O rotational temperatures, of the order of 100 
K, and the order of magnitude gap in the bolometric luminos¬ 
ity between intermediate- and high-mass protostars does not al¬ 
low us to conclude whether the changes are a smooth func¬ 
tion of luminosity. Still, the similarities in rotational tempera¬ 
tures seen for sources with luminosities spanning 6 orders of 
magnitude and probed at different spatial scales strongly sug¬ 
gest the same excitation mechanism, the UV-irradiated shocks 
associated with j ets and winds for all s ources across the lu¬ 
minosity range ()Kristensen et al.l 120 1 3t IKarska et alJ l2014bl 
iMottram et al,ll201417! 
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Table A.l. Log of PACS observations 


Source 

OBSID 

OD 

Date 

Total time 

(s) 

RA 

m s\ 

DEC 

(O , 

AFGL 490 

1342202582 

454 

2010-08-11 

3882 

3 27 38.40 

+58 47 08.0 


1342191353 

290 

2010-02-28 

2029 

3 27 38.40 

+58 47 08.0 

L1641 S3MMS1 

1342226194 

823 

2011-08-14 

3882 

5 39 55.900 

-7 30 28.00 


1342226195 

823 

2011-08-14 

1987 

5 39 55.900 

-7 30 28.00 

NGC 2071 

1342218760 

703 

2011-04-17 

3882 

5 47 04.40 

+0 21 49.00 


1342218761 

703 

2011-04-17 

1987 

5 47 04.40 

+0 21 49.00 

Vela IRS 17 

1342209407 

551 

2010-11-16 

3882 

8 46 34.70 

-43 54 30.5 


1342211844 

594 

2010-12-28 

1987 

8 46 34.70 

-43 54 30.5 

Vela IRS 19 

1342210189 

552 

2010-11-16 

3882 

8 48 48.50 

-43 32 29.0 


1342210190 

552 

2010-11-16 

1987 

8 48 48.50 

-43 32 29.0 

NGC 7129 FIRS2 

1342186321 

165 

2009-10-26 

3895 

21 43 01.70 

+66 03 23.6 


1342186322 

165 

2009-10-26 

2000 

21 43 01.70 

+66 03 23.6 


Appendix A: Supplementary material 

Table IA.il shows the observing log of PACS observations including observations identifications (OBSID), observation day (OD), 
date of observation, total integration time, and pointed coordinates (RA, DEC). Table IA.21 shows molecular data and observed lines 
fluxes and upper limits for all sources that are analyzed in the paper. Figure IA.1I shows selected spectral regions to illustrate the 
quality of the data. Figure [Aj2] illustrates the patterns of continuum emission at 145 pm and the CO 18-17 line emission at 144 pm 
toward all the sources. 






Table A.2. Molecular data" and observed line fluxes 


Species 

Transition 

Wave. 

(pm) 

Rest Freq. 
(GHz) 

EJk B 

(K) 

Aui 

(s’ 1 ) 

AFGL 490 

L 1641 

NGC 2071 
Fluxes (10 

Vela 17 
~ 20 W cm- 2 ) 

Vela 19 

NGC 7129 FIRS 2 

CO 

14-13 

185.999 

1611.8 

580.5 

2.7C-4) 

11.01 ±0.11 

5.00 ±0.10 

100.13 ±0.69 

14.92 ±0.15 

7.38 ±0.10 

4.41 ± 0.07 

h 2 o 

2 2 i-2i 2 

180.488 

1661.0 

194.1 

3.1(-2) 

1.62 ± 0.07 

0.69 ± 0.06 

12.95 ± 0.09 

2.74 ± 0.06 

0.74 ±0.12 

1.28 ±0.05 

h 2 o 

2j 2 -loi 

179.527 

1669.9 

114.4 

5.6(-2) 

absorption 

1.69 ±0.06 

39.81 ±0.40 

9.65 ± 0.08 

2.88 ± 0.06 

3.85 ± 0.08 

h 2 o 

3o3-2i 2 

174.626 

1716.8 

196.8 

5-1(-2) 

4.36 ± 0.08 

2.25 ± 0.06 

36.97 ± 0.30 

6.58 ± 0.07 

3.01 ±0.11 

3.13 ±0.07 

CO 

16-15 

162.812 

1841.3 

751.7 

4.1(-4) 

8.84 ±0.11 

3.76 ± 0.09 

85.61 ±0.45 

9.48 ±0.18 

6.32 ± 0.08 

3.77 ± 0.07 

CO 

18-17 

144.784 

2070.6 

945.0 

5.7(-4) 

7.74 ±0.12 

3.36 ±0.13 

70.00 ± 0.65 

7.62 ±0.18 

5.84 ± 0.08 

3.51 ±0.13 

h 2 o 

3b-2o 2 

138.528 

2164.1 

204.7 

1.3(-1) 

2.86 ± 0.08 

1.57 ±0.07 

25.60 ±0.14 

3.99 ± 0.06 

1.81 ±0.06 

2.23 ± 0.06 

h 2 o 

4o4‘3l3 

125.354 

2391.6 

319.5 

1.7(-1) 

0.89 ±0.11 

1.80 ±0.05 

15.52 ±0.20 

2.19 ±0.10 

1.48 ±0.11 

1.57 ± 0.08 

CO 

22-21 

118.581 

2528.2 

1397.4 

1.0(-3) 

4.61 ±0.15 

2.37 ±0.11 

49.00 ± 0.26 

3.68 ±0.15 

2.78 ± 0.52 

2.70 ±0.11 

CO 

24-23 

108.763 

2756.4 

1656.5 

1.3C-3) 

3.91 ±0.34 

2.25 ± 0.43 

38.80 ± 0.96 

3.30 ± 0.29 

4.80 ± 0.46 

3.23 ± 0.33 

h 2 o 

2 2 i - 1 io 

108.073 

2774.0 

194.1 

2.6(-l) 

1.47 ±0.17 

3.19 ±0.14 

35.30 ±0.30 

5.83 ±0.16 

3.12 ±0.16 

4.51 ±0.12 

CO 

29-28 

90.163 

3325.0 

2399.8 

2.1(-3) 

<0.1 

1.58 ±0.57 

> 8.17 

< 0.46 

1.04 ±0.71 

1.58 ±0.42 

h 2 o 

3 22 -2h 

89.988 

3331.5 

296.8 

3.5(-l) 

<0.1 

0.99 ±0.57 

10.21 ±0.30 

< 0.46 

<0.20 

1.25 ± 0.42 

CO 

30-29 

87.190 

3438.4 

2564.9 

2.3(-3) 

< 0.15 

1.58 ±0.09 

20.37 ±0.13 

2.19 ±0.19 

2.22 ±0.16 

1.58 ±0.07 

h 2 o 

7,6-707 

84.767 

3536.7 

1013.2 

2.1M) 

<0.1 

< 0.10 

< 1.09 

< 0.13 

<0.16 

< 0.10 

CO 

32-31 

81.806 

3664.7 

2911.2 

2.7(-3) 

<0.1 

0.92 ±0.57 

13.16 ±0.40 

< 0.24 

< 0.21 

0.96 ± 0.29 

CO 

33-32 

79.360 

3777.6 

3092.5 

3.0(-3) 

< 0.13 

0.48 ± 0.07 

12.56 ± 0.24 

< 0.10 

<0.16 

1.14 ±0.09 

h 2 o 

6l5“5 2 4 

78.928 

3798.3 

781.1 

4.6(-l) 

< 0.17 

0.53 ±0.09 

2.71 ± 0.43 

< 0.16 

<0.23 

< 0.16 

CO 

36-35 

72.843 

4115.6 

3668.8 

3.6(-3) 

< 0.15 

< 0.04 

7.74 ± 0.21 

< 0.18 

<0.36 

0.70 ± 0.09 

CO 

38-37 

69.074 

4340.1 

4080.0 

4.1(-3) 

< 0.16 

< 0.11 

6.36 ± 0.35 

< 0.13 

<0.35 

< 0.17 

HtO 

OO 

oo 

o 

'J 

63.324 

4734.3 

1070.7 

1.8(0) 

< 0.27 

0.64 ±0.11 

8.36 ± 0.37 

< 0.15 

<0.23 

0.57 ± 0.07 


Notes. Compiled using the CDMS dMuller et alJl200lL 12005 1 and JPL llPickett et alJ 19981) databases. Values of Einstein A coefficient are written in a form A(B) = Ax 10 B . 
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Fig. A.l. Spectral scans covering selected FbO, CO and OH lines in the intermediate mass protostars from the WISH program. The 
rest wavelength of each line is indicated by dashed lines: blue for H 2 O, red for CO, and light blue for OH. 
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Fig. A.2. PACS spectral maps in the CO 18-17 line at 144 pm and the continuum emission at 145 pm in orange contours corre¬ 
sponding to 30%, 50%, 70%, and 90% of the peak value written in the bottom left corner of each map. Wavelengths in microns are 
translated to the velocity scale on the X-axis using laboratory wavelengths (see Table IA,2b of the species and cover the range from 
-600 to 600 kms~'. The Y-axis shows fluxes in Jy normalized to the spaxel with the brightest line on the map in a range -0.2 to 1.2. 

































































































































































































































































































































































































































































































































































































